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Abstract. Galileon models are a class of effective field theories that have recently 
received much attention. They arise in the decoupling limit of theories of massive grav- 
ity, and in some cases they have been treated in their own right as scalar field theories 
with a specific nonlinearly realized global symmetry (Galilean transformation). It is 
well known that in the presence of a source, these Galileon theories admit superluminal 
propagating solutions, implying that as quantum field theories they must admit a dif- 
ferent notion of causality than standard local Lorentz invariant theories. We show that 
in these theories it is easy to construct closed timelike curves (CTCs) within the naive 
regime of validity of the effective field theory. However, on closer inspection we see that 
the CTCs could never arise since the Galileon inevitably becomes infinitely strongly 
coupled at the onset of the formation of a CTC. This implies an infinite amount of 
backreaction, first on the background for the Galileon field, signaling the break down 
of the effective field theory, and subsequently on the spacetime geometry, forbidding 
the formation of the CTC. Furthermore the background solution required to create 
CTCs becomes unstable with an arbitrarily fast decay time. Thus Galileon theories 
satisfy a direct analogue of Hawking's chronology protection conjecture. 
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1 Introduction 

Mankind has long dreamed of being able to travel back in time, although the discussion 
of the associated causal paradoxes has an almost equally long history. The unification 
of space and time in Einstein's relativity led to the possibility of considering time travel 
as a physical possibility that could be studied mathematically. 

As human history has not been invaded by multitudes of visitors from the future, 
there is strong evidence that travel into the past cannot occur, however the question 
of whether the laws of physics allow the possibility of traveling back in time remains. 
Indeed it is surprising that General Relativity (GR) has perfectly valid solutions which 
allow for Closed Time-like Curves (CTCs). These are geodesic or non-geodesic curves 
which form closed loops in space and time which can be traversed in a time-like manner. 
If a person were to travel along such a path they would arrive back at their initial point 
in coordinate space and time after a (positive) finite amount of proper time had elapsed. 
Just as the ideas of time travel immediately raise the ideas of causal paradoxes, the 
existence of a CTC makes the causal structure of a spacetime impossible to determine. 
There is no Cauchy surface on which initial data can be given which determines the 
future of the spacetime uniquely. 
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In 1991 Hawking argued that physics contains a mechanism which protects itself 
against the formation of CTCs, [1], (see Ref. [2] for a comprehensive review). The 
specific example he used was the two-dimensional Misner geometry however he argued 
that the qualitative properties of this model may hold in general. If a CTC exists 
in spacetime the energy-momentum of the quantum field which describes the particle 
trying to traverse this path becomes so large that its backreaction cannot be neglected, 
this is then expected to destroy the existence of the CTC by modifying the spacetime 
geometry itself. In modern parlance, the effective field theory of gravity inevitably 
breaks down before the onset of formation of the CTC. This was called Chronology 
Protection by Hawking. 

In this paper we explore how chronology protection occurs in the Galileon model 
and massive gravity. One key difference worth emphasizing at the outset is that in 
these models there arises more than one effective metric describing the propagation of 
different species, specifically the effective metric seen by the Galileon, and the usual 
metric of the graviton. GR contains a single metric which imposes causality, and con- 
ventional matter coupled to GR is (sub)luminal, nevertheless it does admit CTCs for 
sufficiently peculiar geometries. In Galileon models and massive gravity, the situation 
is seemingly worse since on top of that they allow for fluctuations in the Galileon field 
which are superluminal with respect to the usual GR metric. Thus we anticipate the 
possible existence of new types of CTCs, not allowed in GR, which are tied to the 
existence of this two metric structure, i.e. arise as a direct consequence of this super- 
luminality. In fact we shall show that it is easy to find a situation whereby both metrics 
are well-defined in the sense that particles coupled exclusively to a single metric will 
never be able to form CTCs. However, allowing for the two species of particles to 
interact will allow for an exchange of information between the two metrics which will 
in turn allow for a generalization of the concept of a CTC to arise. In this generalized 
notion of a CTC - 'bimetric CTCs' - it is sufficient that the curve is timelike only with 
respect to one of the two metrics at any given proper time along the curve. We shall 
nevertheless see that such generalized CTCs are forbidden by the same chronology 
protection physics. 

A second crucial difference between the CTCs considered in this paper and those 
considered in GR is that they arise even in the decoupling limit M p i — > oo in which 
limit the graviton (GR) metric is taken to be Minkowski g^ u = rj^. In Hawking's 
chronology protection it is the backreaction of the stress energy on the GR metric 
which is seen to forbid the CTC. However, in the present case this backreaction can be 
made arbitrarily small by making Mp\ arbitrarily large. Nevertheless we shall see that 
even before the backreaction on the metric can become significant, the effective field 
theory of the Galileon will breakdown because of strong coupling. In other words, the 
breakdown of the Galileon effective field theory (EFT) will always occur before, i.e. 
at lower energy scales than, the equivalent breakdown in GR. This is tied to the fact 
that the characteristic scale of interactions in these theories A is much lower than the 
Planck scale A <C Mpi. Despite this we shall see that there is an entirely analogous 
(but for these reasons not identical) chronology protection mechanism in these models 
to the pure GR case discussed by Hawking. 
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The Galileon [3] was proposed as a generalization of decoupling limit of the Dvali- 
Gabadadze-Porrati (DGP) model [4], the first theory of gravity mediated by an effec- 
tively softly massive spin-2 field without incurring ghost-like pathologies. In the de- 
coupling limit of DGP, [5], where the Planck mass becomes infinite and the mass of the 
graviton vanishes in such a way that the spin-2 and spin-0 components of the massive 
graviton decouple, the resulting effective scalar field theory has two defining properties: 
The theory obeys the Galileon symmetry 7r — > it + v^x^ + c, and the equation of motion 
for the scalar field is second order in derivatives despite higher order derivative terms 
appearing in the Lagrangian. Nicolis et.al. showed that in four dimensions there are 
five Lagrangian operators with these two properties and they dubbed the most general 
scalar field theory of this form the Galileon [3]. 

In its original formulation the Galileon was expected to arise within the effective 
description of a more fundamental theory, similarly as in DGP. However it was soon 
considered within the community as a fundamental scalar degree of freedom in its own 
right, with no reference to another underlying description [6-8]. Only later, was it 
realized that the Galileon generically appears as the helicity-0 mode in the decoupling 
limit of healthy theories of massive gravity both in four [9-1 1] and three dimensions 
[12]. Whether or not the graviton has a mass is another unresolved fundamental 
question, and the presence of a graviton mass could have important implications for 
Cosmology. Observations of the solar system constrain the mass to be less than or 
of the order 10~ 32 eV (this bound may change slightly depending on the two free 
parameters of the theory), but this could still become relevant on Hubble scale distances 
today and modify the evolution of the Universe. It has been successfully shown that 
the construction of a massive graviton proposed in [11] is free of ghosts in the ADM 
formalism, first to fourth order [11] and then in full generality, [13, 14]. These results 
have been confirmed in both the Stiickelberg and helicity-languages, [15, 16]. The 
Galileon also arises as the four dimensional low energy effective field theory of certain 
five dimensional probe brane scenarios [7], where the Galileon scalar describes the 
position of the brane in the fifth dimension. Similarly they arise naturally in the 
decoupling limits of higher dimensional extensions of the DGP model [17, 18]. However 
it is the massive graviton origin of the Galileon which is of most interest for this work. 

The Galileon has a broad and interesting phenomenology, but one potentially 
worrying phenomenon is that fluctuations of the Galileon field can propagate super- 
luminally [19, 20]. Superluminal fluctuations had also previously been shown to be 
present in the DGP model [19, 21]. Indeed superluminal modes are generic to Galileon 
[22] and massive gravity [23-25] constructions. Superluminal fluctuations are often 
considered to be a symptom of a sick theory since if a particle is traveling superlumi- 
nally in one reference frame, then there exists another choice of frame in which the 
particle is traveling backwards in coordinate time, and from this a CTC could possibly 
be constructed. However it is not guaranteed that CTCs can form [26]. In general the 
field that exhibits superluminal fluctuations comes with its own effective metric and 
causal structure which is independent to that felt by photons. Even if the causal cones 
of these fluctuations lie outside the causal cones of photons, the causal structure of the 
spacetime can be protected [26] if there exists one foliation of spacetime into surfaces 
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which can be considered as Cauchy surfaces for both metrics. 

Fluctuations of the Galileon fields do carry their own effective metric and hence 
have an independent causal structure. The superluminality arises when there is a non- 
trivial background configuration for the Galileon ttq(x), for which fluctuations around 
this background configuration feel a metric G^ v which depends on the metric of the 
background spacetime, but also on ttq(x) and its derivatives. For suitable choices of 
7r (x) this can allow for superluminal Galileon fluctuations. 

It was suggested recently that CTCs [27], could also be constructed in the Galileon 
model. It was shown that Galileon fluctuations propagate superluminally when the 
background is a plane wave solution of the Galileon equations of motion. It was then 
hypothesized that arranging for four of these plane waves to form a rectangle would 
allow for the fluctuations to travel on a CTC. In section 4 we will consider a related 
scenario where the plane wave travels on a circle, and discuss the circumstances under 
which a CTC forms, and how the Chronology Protection Conjecture operates in this 
scenario. 

In this paper we consider the formation of CTCs in two different setups. The 
first one is unique to massive gravity, whilst the second one relies on the existence of 
two different background metrics and could arise both in massive gravity and in more 
general Galileon theories. In both cases, we show explicitly that the effective metric 
felt by the graviton helicity-0 mode or the Galileon is identical to a well-known class 
of metrics of GR for which CTC may be created, and an analogue of the Hawking 
Chronology Protection Conjecture is found. In particular we see that when starting 
from healthy initial conditions, one necessarily needs to pass through a region of infi- 
nite strong coupling to produce the required background on top of which CTCs may 
form. Furthermore we emphasize that modes with large enough momentum along the 
direction transverse to the CTC are unstable. Since one cannot prevent the excitation 
of such modes at the quantum level, this implies the existence of an arbitrarily fast 
instability. 

The rest of this paper is organized as follows: We start by reviewing the formation 
of CTCs in GR and Hawking's Chronology Protection Conjecture. We then move to 
the description of the decoupling limit of massive gravity in section 3 and as an explicit 
example we show how a simple gas of dust may be sufficient to produce a background 
solution for which CTCs may form. We show that precisely the same arguments as 
in GR hold in the case of massive gravity, and one can therefore extend Hawking's 
Chronology Protection Conjecture to massive gravity for configurations where the ef- 
fective metric exhibits CTCs even if the real metric does not. We then present another 
class of configurations in section 4 which are applicable for both Galileon and massive 
gravity models. In that case CTCs may form when particles living on two different 
effective metrics interact. However in the regime where CTCs form, the background is 
unstable and decays arbitrarily fast. Furthermore, it is impossible to reach this regime 
without the Galileon becoming infinitively strongly coupled. This implies that the 
CTC could never have formed, at least within the regime of validity of the EFT. We 
finally conclude in section 5 by summarizing these results and providing an outlook. 
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2 Closed Time-Like Curves in General Relativity 



As discussed in the introduction CTCs are known to exist in General Relativity. Here 
we review one example and explicitly demonstrate how the chronology protection oc- 
curs. This example, although by itself somewhat over simplistic, will be closely related 
to our later construction of CTCs for the Galileon and massive gravity. It does however 
exhibit many of the generic features expected of CTCs. 
We consider the cylindrical metric, [2] 



ds 2 = -(dt + a(x)dx) 2 + dx 2 + dy 2 + dz 2 



(2.i; 



where x is a periodic coordinate and we identify x = x + L. The function a(x) is an 
arbitrary periodic function. We can always perform a coordinate transformation to 
replace a(x) with A = L~ x J Q dxa(x) by means of t — > t — f(x) with f'(x) = a(x) — A, 
however it is never possible to remove the constant A because the associated coordinate 
transformation f(x) would not be periodic. 

The null vectors which point along the spatial x direction are 



/ 1 — a\ 
1 





V 2 



fl + a\ 
-1 


V J 



(2.2) 



It is clear that for sufficiently large \a\ the causal cones tips over to include points at 
previous values of coordinate time. Propagating backwards in coordinate time at one 
instant is not enough to create a CTC. However this metric can give rise to CTCs for 
a suitable choice of a(x) or more precisely A. 

For concreteness, let us consider a "right-mover" particle. If it starts from the 
origin at t — 0, it will reach the point x > after a time, 



t(x) = ax — / a(x)dx 



(2.3) 



where a is a constant. This is a null geodesic if a = 1, and timelike if a > 1. The 
geodesic path returns to the same point in space after traversing the interval 16 [0, L] 
at time 



T = aL 



a(x)dx . 



(2.4) 



and this occurs at negative coordinate time if a is such that the following condition is 
satisfied 



1 < a < 



1 



a(x)dx = A 



(2.5) 



This condition can only be satisfied if there exists at least one point in the interval 
x G [0, L] for which a(x) > 1 (this is a necessary but not sufficient condition). For a 
time-like curve, increasing a makes this bound tighter. 
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2.1 Instabilities 



So it is possible for CTCs to exist in such a metric, but the Chronology Protection 
Conjecture teaches us that physics can protect itself from particles traveling around 
these paths. To see that this is exactly what happens for this metric consider the 
action for a massless scalar field on this background with action 

S = Jd 4 x^ [(1 - a 2 )^) 2 + 2ad t j>d x <j> ~ (d^f - (d^f - (d z <f)) 2 } , (2.6) 

and equation of motion 



(1 - a 2 )d 2 <f) - 2ad x d t <f) - a'{x)d t <j> + &J> + d 2 <p + d 2 z <f) = 0. 



(2.7) 



The general solutions can be written as superpositions of eigensolutions of the form 
<f) — 0o exp ^— iEt + zkx + ik.y — %E J a(x)dx^j , (2.8) 

where from now on we use the notation k = {k y , k z } and y = {y, z}. </>o, E, 
constants that satisfy the familiar mass-shell condition 



E 2 = k 2 + k 2 



(2.9) 



with k 2 = \k\ 2 = k 2 + k 2 . Solutions to the equation of motion must be periodic under 
x — > x + L therefore the momentum along the x-direction must be quantized, and so 
is the energy: 

- f L 

2im = n n L — E(n, k) I a(x)dx = n n L — E(n, k)A L , with n G Z. (2.10) 

Jo 

Combining (2.9) and (2.10), gives an equation for the quantized allowed energies E n 

1 2vm 



E,(n, k) 



'I- A 2 



2itn 



[1-A 2 )k 2 



(2.11) 



When the field is quantized, we should split the modes up into positive and negative 
energy solutions. In this case, this corresponds to a different choice of the sign in front 
of the square root. These satisfy J5 g= _i(— n) = — E q=+ i(n). From now on we focus on 
the positive frequency mode only (q = +1) and drop the subscript for simplicity (i.e. 
E = E q=+ i). Putting this together, the correctly normalization quantized fields are 

1 °° f d 2 k ( \ 
0(t,x,y) = ~Yl I T^y { a n,k u n,k( x ^y^) + a l,fe<,fe(^^*)j > ( 2 - 12 ) 



where 



a nrk ,al p ] = (2n?L6 nn ,6W(k-k') 



(2.13) 
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and 



u n j^(x,y,t) = — exp ( —iEt + in n x + ik.y — iE f a(x)dx J . (2.14) 

J2(E - AK n ) \ Jo J 



^2{E - An n ) 
The unfamiliar factor 



(A 2 2 \ 
^- + k 2 (l-A 2 )\ (2.15) 

arises from the Klein-Gordon inner product normalization. 

Clearly when \A\ > 1 the energy becomes complex for fixed n at sufficiently high 
k, meaning both that the solutions become unstable, and that the theory cannot be 
quantized. Interesting had we concentrated on a purely two dimensional model of the 
CTC (i.e. k — 0), it would not have been possible to see this instability. We shall see 
a similar result in the CTCs we find later. Furthermore, the instability gets faster at 
larger k meaning that the instability is always faster than any other dynamical time 
scale in the system. 

In the above we have quantized the fields following the standard canonical pro- 
cedure applied to surfaces of constant t. This is the correct procedure for \A\ < 1 
for which t is a good time coordinate. However for \A\ > 1 surfaces of constant t are 
no-longer spacelike, which is precisely the condition necessary to allow for the forma- 
tion of CTC. Nevertheless, what the above solutions show is that regardless of how we 
choose to quantize in that region, an instability is present since for \A\ > 1 eq. (2.11) 
implies either that E is complex at sufficiently large k or that k y or k z are complex. 
Crucially it is the quantization condition which allows us to conclude this, and in this 
sense this is a purely quantum mechanical effect which would not have been seen by 
simply looking at the classical geodesies. 1 

An obvious criticism of solutions of this type that has been often levied is that if 
|A| > 1 and A is independent of time, then the CTC has existed for all time, this is 
then an example of what has been called the 'garbage in - garbage out' principle, [2]. 
We should only really worry about CTCs if they can be shown to inevitably arise from 
well-defined initial data. However, the key point is that even if we imagine forming 
this CTC from some previously well defined solution, for instance let us imagine A 
evolves in time from the region \A\ < 1 to the region \A\ > 1, no matter how rapidly 
the CTC forms, there are always quantized modes of (p m the full theory which have 
sufficiently large momenta that they vary more rapidly than the function A(t). For 
these modes, it is always possible to perform a WKB approximation, for which the 
approximate WKB energies are as in (2.11) with A now replaced by the appropriate 
function of t. Starting from \A\ < 1, the energy of the solution becomes infinite 
when A 2 = 1. This is exactly the condition that was required for the existence of a 
CTC. At this point the backreaction of the quantum fluctuations of the scalar on the 
background geometry will become infinitely large. This can be seen explicitly by either 



1 Of course the same effect is seen for classical fields, but the interpretation of fields in terms of 
particles is intrinsically quantum. 
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calculating the expectation value of the n-point functions of the scalar or for instance 
by calculating the Casimir energy associated with the compact direction x. Therefore 
we clearly see how the Chronology Protection Conjecture manifests in this scenario: 
If the CTC is assumed to form then particles become unstable (or more precisely the 
notion of particles is ill-defined) and we can no longer trust the EFT description of 
the scalar coupled to gravity. However, before the onset of the formation of the CTC 
the EFT will become strongly coupled due to the infinite energies associated with 
arbitrarily large but finite momenta modes. 

2.2 Strong Coupling and Quantum Backreaction 

The essence of the chronology protection mechanism is that quantum fluctuations 
become arbitrarily large before the onset of the formation of the CTC, preventing 
either the formation, or at least its description within the EFT. To understand how 
the chronology protection works consider the Hamiltonian for the scalar. The conjugate 
momentum for the scalar field is 



a 2 )dt<f> + ad x <f>, (2.16) 



making the Hamiltonian 



n = a _ d ^ + \ {{d x <t>? + (d y <pf + {d z <pf) (2.i7) 

= 1(1 - « 2 )(9 t 0) 2 + \ ((«9,0) 2 + («V) 2 + («9 2 0) 2 ) . (2.18) 

For fixed conjugate momentum p^ whenever a solution crosses a — 1, which we recall 
was a necessary condition for the existence of a CTC, the Hamiltonian diverges unless 
simultaneously p^ = d x (f> at that point. Such a condition would remove one of the 
phase space degrees of freedom, and is invariably inconsistent at the quantum level as 
we shall now show. 

First let us compute the quantum expectation value (O|0 2 (x)|O). As usual, this 
is infinite, however it contains an L dependent contribution which is finite due to the 
fact that the divergent terms can all be removed with local counterterms which are L 
independent. Plugging in the quantum modes, and computing in the standard vacuum 
that satisfies a n g|0) =0 we find 

( O |0 2 (x)|O) = - Y f d \ - 1 . (2.19) 

To compute this let us use zeta function regularization (including the parameters \i 
and s — > 0) and replace with 

(O|0 2 (:r)|O> = lim - J (^2^ ((27rn/L)2 + k2{1 ~ A ^y l ' 2 ~ S ■ ( 2 ' 2 °) 
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The integral over k may easily be performed by going to polar coordinates in the y — z 
plane to give 

l-2s 



(O|0 2 (x)|O) = lim 



1 



2n(l - 2s) L(l - A 2 ) 



_1 V 1 s 

I -A 2 ) ^ 2^ 



27m 



L 



(2.21) 



Now in the limit s — > and using the fact that the Riemann zeta function satisfies 
C(-l) = -1/12 so that Er=-oo M = 2 C(~1) = -1/6 we obtain 



(O|0 2 (x)|O) 



1 



12 L 2 (l — A 2 ) 



(2.22) 



Having removed the quadratic divergence, we should really think of this as a finite 
contribution to this correlation function 2 . Its importance is that it will show up in 
interactions, for instance, in the mean field approximation we would replace a A</> 4 
interaction with a term of the form |A(O|0 2 (x)|O)0 2 . More generally it contributes to 
the Feynman propagator any so this L dependence shows up in any perturbative cal- 
culation. The crucial point is that as \A\ approaches unity from below, the two point 
correlation function diverges. Since this term is L dependent it cannot be absorbed by 
a local counterterm. The fact that the correlation function diverges even after renor- 
malization indicates that any self interactions will become infinitely large indicating 
strong coupling of the field. However, even if this field is assumed to be exactly free 
from self-interactions, we cannot switch off the interactions of it with gravity. Thus the 
backreaction of it on the geometry is similarly divergent. To see this let us calculate 
the expectation value of the Hamiltonian. This serves as an indicator of the magnitude 
of the backreaction of the scalar on the spacetime geometry. We have 



(0|K|0) 



d 2 k 



1 



(2tt) 2 A(E - Ak u ) 



'1-A 2 )E 2 + 



2nn 



(2.23) 



This expression is of course again infinite before renormalization, but contains a finite 
L-dependent contribution which accounts for the Casimir effect in this geometry. By 
performing a zeta function regularization, we may replace this expression with 



s 00 /* 




(2.24) 



As before, first performing the integral over k gives 

2 A 



(0\H\0) = lim 



s^o 8ttL(1 - A 2 ) 2 \3-2s 



s — 



2irn 



L 



3~2s 



(2.25) 



2 In the absence of gravity we could be easily persuaded that this contribution could also be absorbed 
via a renormalization. However, L is physically a nonlocal function of the metric, and just as in the 
usual Casimir effect, it has a physical consequence due to the fact that L can change dynamically. 
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which in the limit using 



n 



3 



2C(-3) 



1/60 gives 



— oo 



< 0|M|0) = S5(M 



L 4 (l — A 2 ) 2 ' 



1 



(2.26) 



As in the case of the two-point function, this energy density becomes infinite at the 
onset of formation of the CTC. Once the associated curvature R ~ (0|'H|0}/Mp 1 be- 
comes of order Mpj the effective field theory of gravity breaks down. This clearly occurs 
when L\J\ — A 2 is of order Mpi. Thus in practice this amounts to saying that when the 
proper length of the closed loop is Planckian, we can no longer trust the background 
solution that gives rise to the CTC. 

The metric (2.1) is, after an appropriate coordinate redefinition, precisely the 
same as that presented in Ref. [2] in the toy model of section 19.3.2, with £ there = 
t + J a(x)dx, such that the identification (t, 0) = (t,L) here corresponds to the iden- 



tification (t there ,0) = (t thcrc + LA, L) there. As defined in eq. (19.9) of [2], the 'once- 
through-the-wormhole' interval is then 



This interval measures the proper length around the closed curve. When s 2 > the 
curve is spacelike, when s 2 < is is timelike. We thus see that the two-point function 
of the scalar fluctuations scale as 1/s 2 and the energy density scales as 1/s 4 . 

As discussed in [2] the essence of the chronology protection is that generically 
CTCs, which are assumed to arise at some finite time, arise from a previously closed 
spacelike curve becoming timelike. Continuity guarantees that this can only occur if 
at the transition point the curve is null, in which case the proper distance around it is 
zero. This is the distance captured by s = L\J\ — A 2 . Quantum mechanics, embodied 
in the above Bohr- Sommerf eld type quantization condition guarantees that there are 
quantum contributions to correlation functions and energy densities which scale as pos- 
itive powers of 1/s, at the scale anticipated by the uncertainty principle. At the onset 
of the formation of the CTC we have s = and so all these non-local contributions to 
correlation functions diverge. In practice this means there must be some finite positive 
s at which point the effective field theory in question has broken down. In the present 
case this is the effective field theory of gravity, and the fact that it is breaking down is 
often stated as a breakdown of the semi-classical expansion. 

In short, chronology protection will occur whenever the following two criterion 
are met: 

• The formation of the CTC is associated with the proper distance s around a 
closed loop transitioning from spacelike to timelike. 

• Quantum effects, associated with fluctuations around the loop scale as positive 
powers of 1/s. 



s = 



LVl - A 2 



(2.27) 
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In addition, if the CTC is assumed to form, fields become unstable at an arbitrarily 
fast rate, at least in dimensions d > 2. To prove the chronology protection conjecture 
we would have to prove that the previous two conditions are always met. Unfortunately 
this is a difficult task and so the conjecture has so far been mainly confronted on a case 
by case basis. Nevertheless, in the cases where the CTCs are at least approximately 
geodesies, a condition which does not have to hold, then we always expect an analogue 
of the Bohr-Sommerfeld quantization rule to apply, which appears to be sufficient to 
guarantee the conjecture. 

2.3 Localized CTCs 

The previous example, is clearly very special since it was necessary to assume that 
the topology of the Universe was compact in the x direction, and the CTCs were 
associated with traveling around this compact direction. A more physically interesting 
example is a CTC which is localized in space in an otherwise asymptotically well- 
defined spacetime. This is easy to achieve, the previous example is a natural limit of 
a well-known model with localized CTCs, namely the spinning cosmic string metric: 

ds 2 = -{dt + J&ef + dr 2 + (1 - 5f r 2 d6 2 + dz 2 , (2.28) 

where 2ir5 is the deficit angle, 9 = [0, 27r] and J is proportional to the angular momen- 
tum of the string. In the region of fixed r this metric looks formally the same as the 
previous example. However, now the r dependence allows for the CTCs to be spatially 
localized. A straightforward calculation shows that curves at fixed r and z are CTCs 
when 

r < r=~5 ■ (2 - 29) 

In practice the string has some finite core width Lw and the above metric is only 
valid for r > Lw and CTCs can only be said to exist if J > Lw(l — 5). As in 
the previous case, we can imagine forming the CTC at finite time, by adiabatically 
changing the value of S or Lw or J so that there is a transition from the regime 
for which J < Lw(l — $) for which CTCs are expected to be absent, to the case 
J > Lw{l — 5) when they will be sure to exist at distances r < j^. In this sense it 
should be possible to create CTCs which are localized in space and in turn have not 
existed for all times (i.e. are not eternal time-machines). To understand if chronology 
protection is applicable in this case it would be necessary to quantize the fluctuations 
about the string taking into account what happens inside the core of the string since 
that is crucial to provide the transition barrier J > Lw(l — $)■ This is beyond the 
scope of the present work. There exist several other variations on this theme in the 
literature, notably the Gott and Grant solutions, see [2] for an extensive review. 

3 Galileons as Helicity-0 Modes of Massive Spin-2 Fields 
3.1 Decoupling limit of Ghost-free Massive Gravity 

Although the idea of the Galileon arose out of looking at the decoupling limits of higher 
dimensional braneworld models such as the DGP model [4, 5] and Cascading Gravity 
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[17], it is in the context of four dimensional massive gravity models where the full 
connection between Galileons and infrared modified gravity has been easiest to explore 
[9-11]. The first fully consistent theory of massive gravity in four dimensions was 
proposed in [11]. The degrees of freedom of a massive spin two field can be arranged 
into a massless helicity-2 field, two helicity-1 modes and a helicity-0 field. To see 
this explicitly, one can include four Stiickelberg fields <f) a , and consider non-derivative 
interactions for the tensor = g^ u — ri^d^d^ . For a graviton of mass m, the 
Lagrangian is then 

Mp! , / „ m 2 



£ = ~^V=9 ( R - -jU{g, H) ) . (3.1; 



Then defining )C%(g,H) = 5£ — y/5v — Hv the most general potential U that has no 
ghosts is, [11] 

U(g, H) = -4(W 2 + a 3 U 3 + oM), (3.2) 
where the a n are free parameters, and 

U 2 = [/C] 2 -[/C 2 ], (3.3) 
U 3 = [/C] 3 -3[/C][/C 2 ] + 2[/C 3 ], (3.4) 
U 4 = [/C] 4 - 6[/C 2 ][/C] 2 + 8[/C 3 ][/C] + 3[/C 2 ] 2 - 6[/C 4 ] , (3.5) 

where [. . .] represents the trace of a tensor with respect to the metric g^ u . The absence 
of ghost for this theory has been shown in the decoupling limit in [9-11], fully non- 
linearly beyond the decoupling limit in [13, 14], as well as in the Stiickelberg and 
helicity languages in [15, 16]. 

In the decoupling limit, where the Planck mass is taken to infinity Mpi — > oo and 
the graviton mass tends to zero m — > while the strong coupling scale A 3 = M-p\m 2 
stays fixed, the helicity-1 mode decouples and can consistently be set to zero, and the 
scalar mode takes on a Galileon form, [10]. The tensor H^ u is then of the form 

h,u = ^ (v + |a n /- " ^^ n ^ n ^) » ( 3 - 6 ) 

with Il^j, = d^d v ii. 7i is the canonically normalized helicity-0 component of the massive 
graviton and h^ u = Mp\(g^ u — r]^). In the decoupling limit the Lagrangian reduces to 

^ = ~h^h a , + £ ^riW + ^r^T^ (3.7) 

n=l P1 

where £^ is the Lichnerowicz operator. External sources T^ u scale in this decoupling 
limit in such a way that the quantity T^ u = T^/M-pi remains finite. The coefficients 
ci n are related to the constants a n , with ai = 1 and 

X« = Drrg^ - U, u (3.8) 

= K» - DTrn^ - i([n 2 ] - [n] 2 )^ (3.9) 
x$ = ml - 6[u]ul + 3([n] 2 - [n 2 ])n„ - ^([n] 3 - 3[n 2 ][n] + 2[n 3 ]) ,(3.10) 
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where On = d^d^ir. 

For simplicity we will study a massive gravity theory with a 3 = 0, and then absorb 
the value of a 2 into the scale A and set a 2 = 1. The fields can then be redefined to 
diagonalize the kinetic terms of the scalar and spin two degrees of freedom giving 

c = -\hre$hae + - ^{dnfu-n + ^(<9vr) 2 ([n 2 ] - [n] 2 ) 

+7rf - — dyfidv-nT*" . (3.11) 

The equation of motion for 7r is then 
SDvr - A ([ii^] _ [n]2) + _1_ ([n]3 _ 3[n][n2] + 2[n3]) = _f + ^d^irfr. (3.12) 

As already pointed out in Refs. [10, 28, 29], the additional coupling to matter of 
the form d^d^T^ can play a crucial role and distinguishes models of this type 
from generalized Brans-Dicke models where the Brans-Dicke scalar only couples to the 
trace of the stress energy. In what follows we consider a specific configuration which 
allows not only for the superluminal propagation of 7r fluctuations but as well for the 
generation of CTCs. 

3.2 Closed Time-Like Curves 

Similarly to the pure GR case, we consider a cylindrical background, where the coor- 
dinate x is periodic so that x — > x + L, however unlike the previous section, we can 
here focus on a completely flat geometry g^ u = r]^. The appearance of CTCs is then 
due to the fact that tt fluctuations live on a different effective metric G^ v as we shall 
see in what follows. To create a non-trivial background for the field n, we consider a 
perfectly innocent fluid of (pressureless) dust with energy density p 

^wrl aV ' (3 - 13) 

where a is an arbitrary dimensionless constant. This dust will source the metric 
perturbation, however the metric can still be consistently treated as flat at distance 
scales well below the curvature scale which for this solution is set by 1/m where m is the 
graviton mass (assuming a is of order unity). Thus as long as the size of the compact 
direction L is significantly smaller than 1/m, i.e. mL > 1 we may consistently ignore 
the effect of the metric. 

On the other hand, this source creates a non-trivial configuration ttq(x, t) for the 
helicity-0 mode obtained by solving (3.12) 

A 3 

TT Q (t, x) = (a 2 x 2 + 4atx) . (3.14) 

8 

A point of concern is that the solution ttq does not appear periodic in x. However in 
massive gravity tt is not a fundamental object, this role is instead played by d^dvir which 
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is periodic in the coordinate x for this solution. Stated differently, under a coordinate 
shift x — > x + L the field transforms as tx — > ir — ^(«x + It) — which is simply 
a Galileon transformation. In massive gravity, the system is exactly invariant under 
this transformation and hence the physics is unchanged. This solution is unbounded 
at infinity but again this is not a physical problem since only d^d^ir, which is constant, 
is meaningful. Furthermore since we only consider this configuration in a small enough 
region of spacetime x 2 ,t 2 -C Mpi/A 3 = m~ 2 the background for ir is under control in 
this region. 

We now consider fluctuations around this configuration, it = tt + 0, which have 
the following Lagrangian, 

L = -\{G- x Td v ^, (3.15) 
with the effective metric G^ v satisfying, 

{G -iyu = ^_jL {K »»_ Kr]n+ A_ - KK^ - l - ([K 2 ] - [K] 2 ) fA+^T^ 

(3.16) 

where K^ v = d^d^Q. Note that somewhat uniquely to massive gravity models, the 
effective metric is explicitly a function of the local stress energy as a consequence of 
the dudvirT' 11 ' in the action. 

For the present background configuration the inverse effective metric felt by the 
fluctuations is 



/-l + a 2 -a 0\ 
-a 10 
10 
1/ 



(3.17) 



so the effective metric is identical to that discussed in (2.1) in the case of GR. 
Therefore, completely analogously to section 2, CTCs can exist in the metric felt 
by the fluctuations of the helicity-zero mode. Unlike the situation in GR however, 
this happens whilst the true spacetime metric remains flat. Nevertheless, the same 
argument goes through and as discussed in section 2, to construct such a curve requires 
that \A\ = \1/L j adx\ > 1, which in turn implies that the solution itself is unstable, 
with arbitrarily large instability scale k. 

Furthermore, if one were to start with a healthy configuration for which a 2 < 1 
initially then no CTCs would be present in the initial setup. We could then raise 
the question of whether one could adiabatically change the settings so as to reach 
a regime where CTCs could be produced. For this to happen, we need to increase 
the energy density p adiabatically, or in other words make a a slowly time-varying 
parameter, and include a flux of energy for instance T i(x) such as to locally form a 
lump of matter with a locally increasing in time. CTCs can then arise when a crosses 
the threshold a = 1, but as we have seen in the previous section, such a background 
configuration cannot be constructed without the renormalized two-point function of 
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fluctuations around the background diverging. Thus if we are to start from a healthy 
configuration, one cannot produce a background that allows for CTCs without going 
through a infinitively strongly coupled regime in the process. Therefore we conclude 
that it would not be possible to construct such CTCs in a causal manner. 

3.3 Galileons as Goldstone/Stiickelberg fields versus Fundamental Fields 

In the next section, we move onto more generic Galileon models, where we consider tt 
as a fundamental degree of freedom. The key distinction in this description is that the 
shift and Galileon symmetry are then accidental rather than fundamental symmetries 
of a Stiickelberg field/Goldstone mode as was the case in the previous section. So far 
the fact that only the quantity <9 M <9„7r was physical implied that only that quantity 
had to satisfy the periodicity condition. However when considering a setup where the 
Galileon is treated as its own fundamental scalar field, this no longer holds and one 
needs to impose the periodicity condition on tx as well. The previous configuration can 
therefore not be applied for a Galileon when it does not necessarily play the role of a 
Goldstone boson as in Massive Gravity, and as we shall see below one needs to work 
slightly harder to find a configuration that allows CTCs. 

4 Galileons as Fundamental Fields 

On a flat spacetime the Galileon [3] has the following Lagrangian 

C = -\(d*) 2 - ^ U < d *) 2 + ^(vr) + ^£ 5 (7r) + ttT . (4.1) 
The terms £4(71") and £5(71") are given by 

£ 4 (vr) = {dvf ([II] 2 - [IT 2 ]) , (4.2) 
A(vr) = (<9vr) 2 ([IT] 3 - 3[n][n 2 ] + 2[II 3 ]) . (4.3) 

The c n are arbitrary dimensionless coefficients which are expected to be of order one. 
The purely scalar part of the Lagrangian respects the Galileon symmetry, 71 — > 71 + c + 
b^x^ and is defined up to total derivative terms which are irrelevant in flat space. 

The Galileon operators are non-renormalisable, and so the Galileon model should 
be treated as an effective field theory, valid up to some scale A. However the specific 
operators C n are not renormalised by loop corrections [3, 5, 30]. The classical contri- 
bution of the (n+ l)-th Galileon operator is suppressed compared to the n-th operator 
by a scale which can be written schematically as 

(a A^ 

^x = -^. (4.4) 

This is the dimensionless measure of when the non-linearity of the Galileon kinetic 
terms becomes relevant for classical calculations. There is another dimensionless pa- 
rameter in the theory, which schematically is 

d 2 

«Q = ( 4 - 5 ) 
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which is the parameter suppressing quantum loop corrections [3, 30, 31]. The two 
parameters ac\ and «q are distinct, and this allows for solutions for which the non- 
linearities of the Galileon become important classically ac\ ~ 1, whilst quantum cor- 
rections are still well under control aQ C 1 [30, 31]. Thus non-linear terms can be 
important classically without the breakdown of the effective field theory. 

The equation of motion for the scalar field obtained from the Lagrangian (4.1) is 

[n] + ^ ([n] 2 - [n 2 ]) - ^ ([n] 3 - 3[n][n 2 ] + 2[n 3 ]) (4.6) 
- ^ ([it] 4 - 6[n] 2 [n 2 ] + 8[n][n 3 ] + 3[n 2 ] 2 - 6[n 4 ]) = -t . 

For the solutions we consider the higher order Galileon interactions will all vanish and 
so it is sufficient to study the simplest version of the Galileon theory with C4 = C5 = 0. 
We can therefore absorb the coefficient C3 into A, so as to set C3 = 1. Assuming once 
again a flat cylindrical background metric, where the x coordinate is periodic under 
the shift x — > x + L (with L positive), and some background configuration vr (x), 
fluctuations about this background are described by the Lagrangian 

= -\ {G'T 'OMA - i^P^? ■ ( 4 - 7 ) 

Here again, the causal structure that determines the propagation of fluctuations about 
this background is set not by the background metric g^ u but by the tensor 

where as previously, K^ v = d^d^Q. We will see that under the right circumstances 
this metric G^ v combined with the real spacetime metric g^ allows for the formation 
of CTCs. 

4.1 Galileon Fluctuations 

If the Galileon is treated as an effective field theory with no knowledge of its UV 
completion, then the solution for ir which gives rise to CTCs in the massive gravity 
scenario (3.14) is not acceptable because n is not periodic in the coordinate x. Imposing 
that solutions for ir must be periodic in x makes it more difficult to construct a CTC 
in the Galileon model. The metric controlling the propagation of Galileon fluctuations 
in the simplest scenario with C4 = C5 = was given in Equation (4.8). Unfortunately 
there is no solution for periodic 7r for which G^ u matches the periodic metric we 
discussed before in section 2, so we have to try a little harder to construct a CTC. 

In Ref. [27] it was suggested that a CTC could be created in a Galileon theory 
on a background of plane wave solutions. As plane wave solutions satisfy Dtto = 0, 
this simplifies the form of G^ v . Solutions to the wave equation on a flat Cartesian 
background can have the form 

A 3 

n = Y f(x + t), (4.9) 
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for an arbitrary periodic function f(y) = f(y + L), which is such that /" is dimen- 
sionless. The metric felt by fluctuations around this background solution is then 
given by 



G^dx^dx" = -dt 2 + dx 2 + f" (dt + dxf + dy 2 + dz 2 
and the (t, x) null vectors corresponding to this metric are 



(4.10) 



/ 1 \ 
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f" 

J 



(4.11) 



so that fluctuations are superluminal when f"{x + t) < 0. 

We will build on the principles of the proposal of Ref. [27] to construct a CTC 
for the Galileon. We consider a plane wave that travels in a loop along a periodic 
coordinate. Taking the x coordinate to be periodic we must ensure that the function 
f(x + t) is periodic under the shift x — > x + L. Just as the cylindrical CTC discussed 
in section 2 may be viewed as a special limit of a more localized CTC such as in 
spinning cosmic string, we can view the present example as an approximation to a 
Galileon plane wave traveling on say the trapped photon surface outside a black hole, 
or some equivalent construction where the periodic direction x is replaced by an angular 
direction (the black hole example is not an ideal set up as in a Schwarzschild spacetime 
closed orbits for relativistic particles are unstable). 

4.2 Relevance of Two Metrics 

It is worth pointing out that even though the metric (4.10) allows for superluminal 
modes, the metric itself is still a compactified version of Minkowski spacetime. Working 
in terms of u = t + x and v = t — x, the effective metric is simply 



G^dx^dx" = -dudv + f"(u)du 2 + dy 2 + dz 2 
We can then make the change of coordinates, 

V = v- f{u) , 

such that the effective metric is clearly that of flat cylinder 



G^dx^dx" = -dudV + dy 2 + dz 2 . 



(4.12) 



(4.13) 



(4.14) 



In the change of coordinates (4.13), notice that since f'(u) is periodic, V satisfies the 
same identification as v, V = V — L along with u = u + L. If we were for instance in a 
different situation where /" had a zero mode, say /" = a =constant, then the change 
of coordinates v — > V — v — au would have required us to identify V with V — (l + a)L. 
The resulting flat spacetime would then have been cylindrical in a spacetime direction 
rather than just a purely spacelike direction. 
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The fact that the resulting effective metric (4.14) is simply that of a flat cylinder, 
tells us that by means of this metric alone, one would never be able to produce a 
CTC. However this is a different situation to GR, because (4.10) is not the metric seen 
by all the degrees of freedom living in that geometry, but only by Galileons. Other 
fields which are coupled in a standard way to gravity, will see the flat background 
metric g^ v = rj^. Even though one could never produce a CTC if information was 
propagating on the causal light cone of either g^ u or G^ u alone, one can create CTCs 
by combining both metrics. To achieve this it is sufficient to propagate information by 
means of Galileon fluctuations along the null direction provided by V2 in (4.11) for a 
while, before transferring the information to another particle that propagates on null 
or timelike geodesies with respect to g^ v . It is these 'bimetric CTCs' that we shall be 
concerned with in the following. 

Does a Common Cauchy Surface Exist? 

We expect CTCs to form if it is not possible to foliate in the geometries described 
by and g^ v = rj^ via a common family of Cauchy surfaces. To check if this is 
possible, suppose we perform the following coordinate transformation 



where H is an arbitrary function which is periodic under x — > x + L. We would like 
to use the surfaces r— constant as Cauchy surfaces. For these to describe a suitable 
Cauchy surface they must be space-like with respect to both the background metric 
T)^ and the metric G^ u felt by Galileon fluctuations. 

For simplicity we may look at the induced geometry on surfaces of constant y and 
z. A necessary condition for the surface to be spacelike with respect to r]^ v is if 



t = H(x,y,z,r) , 



(4.15) 



rj pLl/ dx^dx u \t =H (x,y,z,T)-Ay=Az=AT=o = (1 - {d x H) 2 )dx 2 > 



(4.16) 




(4.17) 



G fiu dx fl dx 1 ' ' \ t=H(x,y,z,T);dy=dz=dT=0 > 



(4.18) 



which implies 



(1 - (d x H) 2 + f"(H + x)(l + (d x H) 2 ) dx 2 > , 



(4.19) 



or equivalently, 



(1 + d x H) (l - d x H + f\H + x)(l + d x H)} > 



(4.20) 



Combining equations (4.17) and (4.20) we find 



Since H (x, y, z, r) is periodic in x, there must exist some xq for any time r at which 
d x H(x, y, z, t)\ x=xq = 0. Thus a common Cauchy surface can only exist if, at this point 

l + f"(x o + H(x o ,y,z,T))>0. (4.22) 

Since this relation must hold for all r, a Cauchy surface common to both metrics exists 
only if f"(u) > —1 for all u. Conversely if there exists a u for which f"(u) < —1 then 
there is no possible choice of surface which is spatial with respect to both metrics for 
all times r. As we will see in what follows, CTCs may only form if there exists a point 
u for which f"(u) < — 1 which implies that CTCs may only form when both metrics 
share no common Cauchy surfaces as expected. 

4.3 Bimetric CTCs 

To see explicitly how the CTC forms let us first find the evolution of Galileon fluctu- 
ations. Geodesies moving in the +x direction in the metric obey 

t = x+ f f"(u) du (4.23) 
Jo 

= x + f(t + x) , (4.24) 

where without loss of generality we have chosen to set /'(0) = 0. If T is the time taken 
to travel from x = to x = L, along the "right-mover" null geodesic, then 

T = L + f(T) . (4.25) 

It is therefore clear that it is not possible to construct a curve that returns to the 
point x = L within an elapsed time T = 0. However one can choose the background 
configuration appropriately (i.e. by choosing the periodic function f(u)) such that 
there can be solutions to (4.25) with negative T. Such solutions are not yet CTCs, 
but can be used as the basis of constructing bimetric CTCs in which information 
propagates on both metrics. 

One may consider the configuration where information is carried by Galileon fluc- 
tuations starting from the origin 0(x — 0, t — 0) and going to a point A with spacetime 
coordinates x = L = and t = —\T\, with \T\ < L satisfying 

f\-\T\) = -(L + \T\)<-L. (4.26) 

For this condition to be satisfied, there must be at least one point in the interval 
< u < L — \T\ where f"(u) < — 1. However since / is periodic, there must also be at 
least one point in the interval — \T\ < u < where f"(u) > 1, and we can then infer 
that if one were to continue propagating information with the effective metric G^ u , the 
origin point O would not lie within the future light cone of the point A. However if one 
considered interactions at the point A such that from that point on, information was 
instead carried by a standard massive particle propagating in the metric g^ u = i]^, 
this particle could then remain stationary in space but travel forwards in time to reach 
x = at time t — 0. Then we have created a CTC for information, even if no single 
particle traverses the entire curve. 
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4.4 Stability 



We have established that the background plane wave solution will admit bimetric CTCs 
- CTCs of information - provided we allow for interactions between the Galileon and 
another particle. These arise whenever there exists some region for which f"(u) < 
— 1. Let us now analyze the stability of this configuration. We shall see that these 
configurations are unstable, and that the backreaction becomes arbitrarily large in 
precisely the same sense as in the pure GR example discussed in section 2. 

The equation of motion for Galileon fluctuations about the background it = 



IS 



- 4d v d v <f> - 4/»«9> + d 2 y <f> + d 2 z c 
whose solutions can be represented as superpositions of 







[u,v,y,z) 



e i[k v y+k z z-K„v]g/\ 



where 0o is constant and the function g(u) has the form 

~ k 2 



g{u) = exp 



4k, 



■U - K n f{u) 



(4.27) 



(4.28) 



(4.29) 



As before k 2 = k 2 + k 2 is the total momentum along the transverse directions. Peri- 
odicity of is imposed by the following condition 



27m 



for nGZ 



(4.30) 



the periodicity of / means that this condition is independent of the form of /. The 
Galileon field configuration is therefore apparently stable even when reaching the 
threshold /" < —1. This should actually come as no surprise since the effective metric 
G^y is nothing but Minkowski space in a different system of coordinates. Thus we do 
not expect any instability for fluctuations living on top of this geometry if these fluctu- 
ations are decoupled from the other metric. However as we have seen, the metric G^ 
alone is not sufficient to produce CTCs. To build a closed time-like curve it was nec- 
essary for the Galileon to interact with a massive species which propagates according 
to the causal structure of the other background Minkowski metric r]^. We proceed to 
analyze whether the coupled system is stable when the background Galileon solution 
allows for the formation of closed time-like curves. 



Adding a Coupling to Another Field 

In what follows, we therefore consider two different scalar fields. The first one 
corresponds to the Galileon fluctuation n = 7r (x + t) + 0(x^) where ir is given by 
(4.9), so that lives in the effective metric G^. The second field x nas a mass m x 
and lives on the background Minkowski metric g^ v = t]^. If these two fields interact, 
then generically when expanding around the background configuration it = ttq an 
effective coupling will be generated already at quadratic order. To give an example, 
suppose that the original fields interact via the rather natural Galileon invariant term 
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-~fsX ((dfj,d u 7r) 2 — (On) 2 ). Expanding around the background n = 7r + </>(x M ) we will 
generate already in the free theory a term of the form 

£int = -x/*0, (4.31) 

where in this case the operator = (d^d^nod^dvcj) — D7i D(j)). This operator 
plays the role of an effective kinetic coupling between two fields and more generally, 
to preserve the Galileon symmetry it must contain at least two derivatives acting on 
cf). For the following argument the precise form of this coupling will not matter, only 
that some coupling exists already at quadratic order. Note that couplings like this 
will also be generated by quantum corrections even if not already present at tree level. 
We intentionally consider a quadratic coupling such that both fields interact already 
within linearized perturbed theory. If we were to consider only a higher order coupling, 
one would need to work to higher order to see the interaction arising as well as the 
associated instability. 

The real contribution of the interaction is to regulate the square root that will 
appear in the expression for one of the field momenta in (4.40) and so to impose a 
definite notion of positive and negative frequency modes. This regularization occurs 
independently of the exact interaction we consider. To see this, we begin with the 
action to quadratic order 

C = -i (G- i y V d^d v j> - ^rdpXdvX - \Kx 2 ~xM> (4.32) 
so that the coupled equations of motion are 

- AdM - Af"d 2 v <p + tf<j> + d 2 z <p = tf X (4-33) 

-Ad u d vX + d 2 yX + d 2 zX - m 2 xX = fup . (4.34) 

where //' is the adjoint operator to fi: J & 4 xxfi<fi — J d 4 x(j)ftx- If t ne background 
is translation invariant in the v,y,z directions then on Fourier transforming, the po- 
tentially derivative coupling encoded by the operator fi is replaced with a function 
/j,(k, k,u). Here again, we use the notation where k is the momentum in the v direc- 
tion, k y and k z are the momenta in the y and z directions and where k = {k y ,k z } and 
y = {y, z}, k 2 = k 2 + k 2 . Again the form of this function will not be crucial to the 
following argument. The equations of motion become 

-An 2 f" + k 2 £(«,*,«) ] 
\XJ \ fi(K,k,u) k 2 + m 2 x 

Although solving these equations exactly is extremely difficult, it is always straight- 
forward to solve them in the WKB regime, i.e. for sufficiently high momenta k u along 
the u direction such that the variation of k u is small, namely 

(4.36) 
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Then the standard WKB solution is a good one 

(j) ~ e i f Kdu-iKV+ik.y (4.37) 

where and \ are slowly varying prefactors. Working within that approximation, the 
system of equations has a solution if 

/_ AK 2 f" + k 2 + AKk u ^ \ 

d6 H /* fc 2 + m*+4^J ^ 

This determines the two possible li-modes 



Ik 2 + m 2 - Af"n 2 + aJ(4f"K 2 + m 2 ) 2 + 4/V 
feW = . (4.40) 

OK 

where cr = ±1. To get a sense of what this implies, consider modes of sufficiently 
large momenta k such that both the interacting and mass terms maybe be neglected, 
/i, m 2 -C A; 2 . Then we obtain 

fe? ) = -~+5«(r-^iri). (4.4i) 

4 k z 

Comparing this with the result we would get in the absence of an interaction we see 
that there is a crucial difference, even when taking the limit \i — > (or the interaction 
coupling g — > 0). The momentum along the u direction now not only depends on the 
periodic function f"{u) but also on its absolute value. In the single field case, the 
choice of positive versus negative frequency modes is ambiguous because of the fact 
that the effective frequency depends on /", and can change sign as we traverse the 
x direction. However as soon as the field </> interacts with \i th e causal structure of 
X is imposed upon <fi. The smallest interaction requires us to us to fix the notion of 
positive and negative frequency mode for <p once and for all independently of the sign 
of /". As a result, the absolute value of /" should then be taken into account which 
significantly changes the behaviour of the modes. At high momenta, the quantity that 
can be identified as the frequency or energy of each the mode is 

£?W = — (/" - a\f"\) . (4.42) 

4 k z 

This should be compared with the form of the energies in the absence of interactions 

^LoLinteracting = * + ~ ~ \* U" ~ °f") ■ (4-43) 

Thus we see that the limit /j, — > of the vacuum for the interacting system, is not 
equivalent to the tensor product of the vacua for the non-interacting systems, 



to 

non— interacting 



lim £ (<T) ^ E 

lim (interacting vacuum)^^ 7^ |0)<a <S> |0) x . (4.44) 
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Quantization Conditions 

The energy for the fluctuations become quantized because of the periodicity conditions 
4>{u + L,v — L) = (f)(u, v), x(u + L,v — L) = x{ u i v )i which determine the quantization 
condition for k: 

[ k^du + n n L = 2nn with neZ. (4.45) 
Jo 

This is equivalent to the statement 

(l " ^ £ ir(«)|d«) kI - - \e = , (4.46) 

and is easy to solve giving 

^=2(1^) (p^ + ^ll^) +kHl-aA)\, (4.47) 
where we have defined the positive constant, 

A=^£\f"(u)\du>0. (4.48) 

The choice of <j = ±1 in the expression for n n is an independent sign choice from 
the one denoted by a. Choosing the correct combinations of signs is important in 
determining the correct split into positive and negative frequency modes in order to 
quantize the system. The choice that corresponds to the positive frequency /energy 
modes is ^ = +1. This equation admits a real solution for both sets of modes (i.e. 
both a — ±1 ) only if the background satisfies the condition 




1 



L 



A= 2LJ l/»l dn<1 - ( 449 ) 

Remarkably this condition is violated as soon as we enter a regime where CTCs may 
form, as is shown explicitly in the Appendix. Although not immediately obvious, this 
result follows from the periodicity of the function f(u). Therefore as soon as we are 
in such a regime where CTCs may form, modes with large enough momentum k have 
complex k's and in turn complex energies. The solution is therefore unstable with an 
arbitrarily small time scale set by the momentum k. This is similar to a ghost-like 
instability. If one were to start with an initial situation where (4.49) is satisfied and 
adiabatically increase f"(u) so as to break that bound, the momentum n would diverge 
independently of the transverse momentum k. This is the indicator of strong coupling, 
as we shall see in section 4.5. 
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Eigenstates 

The coupling of the two modes means that the eigenstates of the system are a mixture 
of both (j) and \ modes. The eigenstates are determined by solving 



. (4.50) 



An 2 f" + k 2 + Ank^ 

H k 2 + m 2 x + Ank [ u 

Using the solutions for k^ we find that the relevant eigenstates of this system are 




;(CT) = / 2k 2 f" + \m\ - a^{2K 2 f" + \m 2 x ) 2 + ^ 



(4.51) 



To illustrate the behaviour of these eigenstates more clearly we treat the interaction 
terms /i, [p as small, then 

kW = / 2k v + x - *i*r + XI \ + c(wjt) _ (4 52) 

The behaviour of the eigenstates is clearly determined by the sign of 2k 2 /" + |^i 2 , if 
we are working at low momentum (k m x ) then this is dominated by the mass term 
which is always positive, and to first order in fi 

,,<+■> =(;), ,<-»^( m ;). (4.53) 
At high momenta, to first order in /i 

„ W= ( 2 « 2 (/"-*l/"D). (4 . 54) 

If the top component of this vector is non-zero then the bottom component can be 
taken to be zero in the limit fi — > 0. However, when the top component is zero, the 
bottom component cannot be neglected and can be normalized to unity. Given this we 
see that when /" > the eigenstates have the same form as at low momenta but when 
/" < 0, which is precisely when we have superluminality, the roles of the eigenstates 
switch. Thus we have: 

• At low momenta, the light state a = — 1 gives the eigenstate which is mostly 0, 
and inherits its superluminality, and the heavy state o = +1 gives the eigenstate 
which is mostly x an d (sub)luminal. 

• At high momenta, whether the eigenstate for a given choice of a is mostly (su- 
perluminal) <fi or mostly ((sub) luminal) x depends on the local sign of /" and 
thus oscillates in time and in space. 
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4.5 Strong Coupling and Quantum Backreaction 

We have established that by means of two metrics, one can produce CTCs only if there 
exists a point for which /" < —1. However this is precisely when the configuration 
is unstable due to the perturbations having complex energies. To further solidify the 
unphysical nature of the CTCs we will now see that the EFT inevitably breaks down 
before the CTC has a chance to form. 

To see this we need to quantize the coupled fluctuations of both fields. The free 
Lagrangian for perturbations is given by 



2' y<T ' 



^ = ^(1 + /") W) 2 - f'd^dA ~ J(l - f')(^0) 2 - -(^0) 2 - -(<9 Z 0) 2 

1 



—^T^X^vX ~ 2 m x* - X A • (4.55) 

In order to compute the conjugate momenta it is necessary to specify the dependence 
of the coupling operator (x on time derivatives. We shall assume for simplicity that 
it is independent of time derivatives so that it does not contribute to the conjugate 
momentum of either field. In this case the conjugate momentum associated to 4> is 

p, = (l + / , W-/ ,/ 0*0, (4.56) 

whereas that for \ is conventional; p x = d t X- Then the Hamiltonian is 

U = 2(1 + /") ^ + ^) 2 +^x+^( 5 ^) 2 +^( 5 ^) 2 + ^(Vx) 2 +^ 2 X 2 + X/ i0 . (4.57) 

From a classical perspective it is clear from this that the Hamiltonian diverges when /" 
tries to reach the threshold /" = —1, which is a necessary condition for constructing 
a CTC, provided we assume that p^ + d x <p ^ 0. However on the other hand, if we 
can arrange p^ + d x <p to vanish when /" hits —1 then no obvious problems arise. To 
determine if there really is a problem we must calculate the quantum expectation value 
of the Hamiltonian. 

To quantize the fields we must split them up into positive and negative energy 
modes. However since the fields are coupled each field will be composed of the two 
independent positive energy solutions corresponding to the two choices of the sign a. 
The light eigenstate (the one that is mostly 4> - and gives pure <j) fluctuations in the 
limit m x — > oo) corresponds to the choice a = — 1 so that we have 



4 _) = K n,7=+l = W77—^ I —f- + \l ( —T ) +kHl + A)\, (4.58) 




and the energies are given by 



2k 2 + m\- Af"K { - )2 - ^(Af"K { n h + m x ) 2 + A^ 



£H = = «(") + *—LJl * ^ . (4.59) 
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On the other hand the massive eigenstate (mostly x) corresponds to the choice a = +1 
so that 




4+) = k^) = ^_^|=^1 + A /(=^] + k»{l-A)\, (4.60) 
and the energies are given by 

= ^l+t } = 4 +) + • (4.61) 



n 



To confirm that the correct choices of sign have been made in (4.58) and (4.60) we can 
take the limit f"(u) — > and \x — > in which case we get the familiar result for the 
energy of a massless field in lightcone coordinates E(~> —>■«;+ k 2 /4n along with that 
for a massive field 4k + (k 2 + m^/An. 

It is apparent from the above expressions that it is the massive eigenstate a = +1 
and not the light one whose energies become complex in the region of the formation 
of the CTC. Thus it is the backreaction from the massive eigenstate which will im- 
ply strong coupling and prevent the formation of the CTC. This result is somewhat 
surprising since it is the light eigenstate which is superluminal at low energies. How- 
ever, the backreaction comes from the high energy behaviour of the modes where both 
eigenstates oscillate in time between superluminal and (sub) luminal fluctuations. 

The Klein-Gordan inner product for the normalization of the mode functions is 
the sum of the two contributions from each field. Only this combined sum is conserved 
in the interacting case. 



(0Y10X) = -i / d 3 x ir ((i + f)d t <j> - fd x <j>) - ((i + f)d t r - f"d x r) 0] - 

J t=constant 

z / d 3 x [ X '*d tX - d tX '*x] ■ (4-62) 

J t=constant 

Light Eigenstate Fluctuations a = — 1 

In the light eigenstate can be expressed as a superposition of modes of the form 

u n% {x,y,t) = ^-»(/s ( -H«,»l)d«)+^ ) (n,t)(u-» ) )+ft.S ) ( 4 g 3 ) 

as well as its complex conjugate, and similarly for x- Unfortunately, although providing 
a complete set, this set of modes is not orthogonal on the constant t hypersurfaces. To 
see this let us consider the interactions \i to be negligible and the low momenta regime 
where the eigenstate is mostly 0. The nontrivial form of the momentum conjugate to 
'Pi P<P = (1 + f")9t<P ~ f"9 x <t> implies that the Klein-Gordan inner product takes the 
form 

(u'\u) = -i [ d 3 x [u» ((1 + f")d t u - f"d x u) - ((1 + f )<V* - f"d x u'*) u] . (4.64) 

J t=constant 
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Substituting in modes with different values of n we obtain 



(n\k'\n,k) = ^-{2n) 2 5^\k' -k) t dx\ (— + —) + (1 + f"(x))( Kn + /w) 

2 JO { L 4 V K n 

x exp ^2i(4 _) - k { ~ ] )x - i dx(E { -\x,n, k) - E ( ~\x,n', k))^j 1,(4.65) 

which in general does not vanish for n ^ n' even though it does vanish for different k. 
Note that the t dependence does drop out because the spatial x integral is around a 
closed loop, the starting point of the integral is irrelevant and so we can replace x + t by 
x within the integrand. This confirms that this is the correct conserved inner product. 
The fact that the modes are not orthogonal is not surprising because the modes we 
have defined were obtained by Fourier transforming in the v direction, and so will be 
orthogonal when using a definition of inner product on a constant u hypersurface but 
not on constant t hypersurfaces. Unfortunately this makes the direct quantization of 
modes using the constant t Cauchy surfaces rather complex (although by no means 
impossible - since the set is complete there is no problem in principle). 

Fortunately it is relatively easy to bypass this problem by means of the following 
trick. Define the new coordinate 

v = v- | u d«Qr(u)+iir(fi)i) • (4.66) 

Without the modulus sign this is the same transformation we used before to show that 
the plane wave geometry is equivalent to Minkowski spacetime Eq. (4.13). Working 
at sufficiently high momenta that we can neglect the mass m x of the second field (or 
indeed simply assuming the second field is massless) then the energies of the light 
eigenstate in the limit \i — > take the form 

£Hai H + 5!±^in. (4 . 67) 

4/4 ' 

In this limit it is easy to see that the mode functions simply become 



n,k 



u_ Ax,y,t) = 0exp y-i (^~jzj^j u ~ iK L )v + ik.^j , (4.68) 



which are precisely the mode functions for a massless particle in Minkowski spacetime 
with metric 

dsjL) = -dudV + dy 2 + dz 2 . (4.69) 

That is because in this limit, eq. (4.69) describes the correct effective metric G^J that 
fluctuations of the light eigenstate see. In other words, the light eigenstate continues 
to behave effectively like a massless field on Minkowski spacetime. The fact that it is 
necessary to perform a non-analytic coordinate transformation to see this Minkowski 
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metric, is capturing the fact that the effective metric is really a combination of the 
metric seen by and that seen by \. Both of these are individually Minkowski, but 
not in the same coordinates. 

The coordinate redefinition in eq. (4.66) that we performed to obtain this effective 
metric is not analytic, but it is once differentiable which is sufficient for what follows. It 
can always be considered as the limit of an analytic transformation. To make this clear 
we undo the coordinate transformation, then the metric that describes the propagation 
of fluctuations of the light eigenstate is 



d^_ } = _dndi; + \{f\u) + \f"(u)\)du 2 + dy 2 + dz 2 . (4.70) 



We can choose to regulate the metric in a manner consistent with the \i — > limit 
under which is was obtained 

dsj_ } = -dudv + -(f"(u) + y/f"(u) 2 + e 2 )du 2 + dy 2 + dz 2 . (4.71) 

This is still a plane-wave geometry, and all curvature invariants vanish for this metric. 
Thus there is no problem taking the limit e — > 0. 

In this effective Minkowski space that describes the fluctuations of <fr, we can 
define a new time coordinate T = |(u + V) and a new space coordinate X = \{u — V). 
The identification which defines the period of the compactified dimension which in the 
original coordinates takes x — > x + L, with t left unchanged, in the new coordinates 
takes u — > u + L simultaneously with 

V -L(l + A). (4.72) 

Equivalently we can identify new time and space coordinates T = (1/2) (u + V) and 
X = (1/2) {u — V) in which we identify 

X X + L + l -LA, (4.73) 
T — )• T — -LA. (4.74) 

The interval s that defines the proper length of the closed curve which becomes the 
CTC described in section 2.2 is in the present case 



s = Ly/1 + A. (4.75) 

As we can see s never passes through zero. Thus (at one-loop order) the quantum 
fluctuations of the light eigenstate do not blow up as we approach the region of the 
formation of the CTC. This is consistent with the fact that the energies of the light 
eigenstates always remain real and finite for finite momenta k and n. 
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Massive Eigenstate Fluctuations a = +1 

In the massive eigenstate x can be expressed as a superposition of modes of the form 

V .(x y J) _ ^ e -'(/S (+) («Afc)d«)+ki +) (ii,k)(«-ij)+!£j (4 76) 

as well as its complex conjugate. In this case, in the limit \l — > 0, at low momenta the 
eigenstate is mostly x an d the Klein-Gordan inner product is standard 

( v '\ v ) = -i [ d 3 x [v'*d t v - d t v'*v] , (4.77) 

J t=constant 

and so the different modes are orthogonal and it is possible to employ a standard quan- 
tization procedure. Nevertheless at high momenta, the eigenstates oscillate between 
(f) and Xi as m the case of the light eigenstate, and so to quantize it easier to use the 
previous trick and work with new coordinates. We define 

v=v- f du Qr(fi) - iir(«)i). (4.78) 

Working, again, at sufficiently high momenta that we can neglect the mass m x of the 
second field then the energies of the x field in the limit \i — > take the form 

EM „ „(+) + k 2 -2K^\r-\f'\) ^ (4 7Q) 

and in this limit the mode functions simply become 

v n,d x ' M=X ex P (i u + iK{ n ] y + ik.ffj , (4.80) 

which are again precisely the mode functions for a particle in Minkowski spacetime 
with metric 

ds 2 {+) = -dudV + dy 2 + dz 2 = -dudv + \{f"{u) - \f"{u)\)du 2 + dy 2 + dz 2 . (4.81) 

In the present case, in the superluminal region f"(u) < we see that the modulus sign 
acts to double the f"(u) contribution and so in the high momenta regime the massive 
eigenstate is superluminal with respect to the background metric, even though it always 
travels along its own effectively Minkowski metric. The CTCs are formed from curves 
which are timelike with respect to one of these metrics, but not necessarily both at 
the same time. They should be viewed as curves built out of segments where each 
segment is defined on one or the other of the metrics. These disjoint CTCs describe 
the effective propagation between the two metrics, and the effective metric in Equation 
(4.81) is the one that describes the propagation of that information. 

Finally we can again define a new time coordinate T = ^(u + V) and a new space 
coordinate X — \{u — V) such that the cylindrical identification x — > x + L becomes 

X ->■ X + L - l -LA (4.82) 

f ->f+-LA. (4.83) 
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The proper length of the loop s from section 2.2 is in this case 



s = LVl-A 



(4.84) 



which could have been derived from (2.27) with L — > L — \LA and A — > jzj • Thus 
we may immediately borrow the result of section 2.2 to say that when s approaches 
zero, the two-point function of the massive eigenstate modes diverges as 1/s 2 and the 
expectation value of the total Hamiltonian diverges as 1/s 4 because of the contributions 
coming from the massive eigenstate. The condition for s to remain positive is 



However as already mentioned below Equation (4.49) and proven in the Appendix, this 
condition is violated whenever there exists at least one point for which /" < —1, and 
this is precisely the condition for CTCs. Thus we can conclude that the expectation 
value of the two-point function and the expectation value of the Hamiltonian of the 
system diverge precisely at the point of formation of the CTC, namely when s = 0. 
This is consistent with the fact that the energies of the massive eigenstates become 
infinite, at finite k and n, at the onset of formation of the CTC. Since, even in the 
absence of gravity, this implies an infinite amount of backreaction on the Galileon field, 
we conclude that the formation of the CTC lies outside of the regime of validity of the 
Galileon effective field theory. 

The surprise result of this analysis is that ultimately it is the two-point function 
and stress-energy of the massive eigenstate (the one which is mostly x & t low momenta 
but is an oscillating mix of \ and at high momenta) that diverges at the onset of 
formation of the CTC. This is because in the interacting theory, although for momen- 
tum below the scale m x only the light eigenstate is superluminal when f"(u) < 0, for 
larger momenta above the scale m x the massive eigenstate inherits the superluminality 
in an oscillatory manner. This behavior could never have been anticipated by looking 
at the two fields in the absence of interactions. 

5 Outlook 

The possibility of creating CTCs within the standard framework of General Relativity 
has been well established for more than half a century, beginning with Kurt Godel in 
1949. Allowing for this possibility would dramatically shake up the most fundamental 
principle of physics. Fortunately within the context of GR, Hawking's Chronology Pro- 
tection Principle is widely believed to prevent any such configurations from appearing. 

It has been well established that is is possible to construct EFTs with Lorentz 
invariant Lagrangians that nevertheless allow for fluctuations whose effective metric 
allows for faster than light propagation. This occurs in Galileon models and in massive 
gravity models. In this paper, we have explored the possibility of constructing CTCs 
for Galileons. Depending on whether the Galileons are seen as their own fundamen- 
tal degrees of freedom, or instead as the Goldstone mode of a massive graviton, the 




(4.85) 
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situations in which CTCs may form may change slightly, however the nature of the 
Chronology Protection Principle is always the same and is generically applicable. In 
massive gravity, CTCs may be constructed on a cylindrical spacetime which are iden- 
tical to solutions found in GR. However, just like in GR, we find that these solutions 
are quantum mechanically unstable, and if one were to start with healthy initial con- 
ditions, one would need to cross an infinitely strongly coupled region where the EFT 
is no longer a valid description to achieve a configuration where CTCs may form. 

More generally, even if the Galileon is considered as a fundamental degree of 
freedom in its own right, one can easily find configurations where no single particle 
may propagate along a CTC, however information could ride between different species 
of particles, which feel different effective metrics, and hence create a CTC. In this 
situation CTCs may only form if particles with different effective metrics interact, a 
situation which is inevitable in a gravitational theory where all species of particles 
interact at least via gravity. Whilst these configurations may appear classically stable 
with respect to one metric, we show that as soon as these interactions are taken 
into account, the result changes dramatically and the configuration becomes unstable 
precisely when CTCs may form. Hence within the regime of validity of the theory at 
hand, CTCs are not realized. 

The arguments presented in the paper, do not represent a complete no-go for the 
construction of CTCs, and neither, as yet, does the Hawking Chronology Protection 
Principle. However within a simple class of configurations, we have established that 
the creation of CTCs is inexorably tied to strong coupling issues, making such config- 
urations unreliable using the description at hand. As with the Hawking Chronology 
Protection Principle, much needs to be done to prove the general validity of this ar- 
gument. For instance we have not explicitly shown that it will apply in the case of 
a spatially localized CTC. Nevertheless, there are strong reasons to believe that this 
new type of chronology protection, applicable in the case of theories with more than 
one effective metric and superluminalities, is as generic as that proposed by Hawking. 
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A Bounds of /" 

In this Appendix, we show explicitly that as soon as the Galileon configuration can 
admit CTCs, then the bound (4.49) is violated and the solution becomes unstable with 
an arbitrarily fast decay rate. 

As we have seen in section 4.3, CTCs may form if the periodic function f(u) is such 
that f'(0) = and there exists a point — \T\ for which f'(—\T\) = —(L + \T\) < —L. 
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We start by denoting by uq the minimum of f'(u). For simplicity we take uq in the 
interval —L < u < 0. From the previous bound, we immediately infer that 

f(u )<-L. (A.l) 

We now denote by u n the points — L < Uq < U\ < ■ ■ ■ < it at < where f"(u n ) = 0. 
Knowing that /' is negative and minimal at uq, i.e. that /" is positive in the interval 
[uq, Ui] and changes sign at each point u n , we infer that 

u n +l 

\f"(u)\du = (-l) n (/'(Vi) - /'(«»)) > . (A.2) 

71 

Knowing this, it is easy to see that the total integral within the interval [uo, 0] is already 
greater than L. Denoting un+i = 0, the integral is 

N N 

/ |/»|d« = \f'^)\&u = - f\u n )) . (A.3) 

•> u ° n=0 n=0 

We now first assume that = 2M is even, so that 

/ |r(u)|d« = -f(u ) + 2 Y,(f>2n-i) ~ f(u 2n )) + /'(0) > L . (A.4) 

^"o n= i 

Similarly, if N = 2M + 1 is odd, 

p0 m 

/ |r(«)|d« = -f (n ) +2^(/'(« 2n _ 1 ) - f(u 2n )) + 2f'(u N ) - f(0) > L,(A.5) 

•^"o n= l 

where we used the fact that /'(ttjv+i) = /'(0) = 0, and from (A.2) we have /'(«2n-i) _ 
f{v>2n) > and if A is even, /'(%) > f'(u N+1 ) = 0. 

We could now repeat precisely the same argument within the interval [0, L — uq]. 
If N = M+ 1 is the total number of zeros of the function /" within a period, f"(u n ) = 
for all n — 1, • • ■ , N, then the integral within a period is 

|f (u)|du = -2/'(«o) +2^(/'(« 2n _ 1 ) - f (u 2 «)) + 2fK) > 2L, (A.6) 

-0 n=l 

where we have used the fact that since u = (L — u ) is a minimum, /" is negative in 
the interval [u N , L - u ], so that J u ~ u ° |/"(u)|dtt = f'{u N ) - f'(L - u ) = f'(u N ) > 0. 

We can therefore infer that if the background configuration is in a regime which 
admits CTCs, that is if there is at least one point along the cylinder where f'iu) < —L, 
then this together with the periodicity condition automatically implies that 

[ L \f"(u)\du>2L. (A.7) 
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